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ABSTRACT 

We present 2.5-5.0^m spectra of 83 nearby (0.002 < 0 < 0.48) and bright {K < 14mag) type-1 
active galactic nuclei (AGNs) taken with the Infrared Camera (IRC) on board AKARI. The 2.5- 
5.0/rm spectral region contains emission lines such as Br/3 (2.63/rm). Bra (4.05/rm), and polycyclic 
aromatic hydrocarbons (PAH; 3.3/xm), which can be used for studying the black hole (BH) masses 
and star formation activities in the host galaxies of AGNs. The spectral region also suffers less dust 
extinction than in the ultra violet (UV) or optical wavelengths, which may provide an unobscured 
view of dusty AGNs. Our sample is selected from bright quasar sur veys of Paloma r -Gree n (PG) and 
SNUQSO, and AGNs with reverberation-mapped BH masses from IPeterson et al.l (120041) . Using 11 
AGNs with reliable detection of Brackett lines, we derive the Brackett-line-based BH mass estimators. 

We also find that the observed Brackett line ratios can be explained with the commonly adopted 
physical conditions of the broad line region (BLR). Moreover, we fit the hot and warm dust components 
of the dust torus by adding photometric data of SDSS, 2MASS, WISE, and ISO to the AKARI 
spectra, finding hot and warm dust temperatures of ^ HOOK and ~ 220K, respectively, rather than 
the commonly cited hot dust temperature of 1500 K. 

Subject headings: galaxies: active - galaxies: nuclei - quasars: emission lines - quasars: general - 
infrared: galaxies 


1. INTRODUCTION 

An active galactic nucleus (AGN) is a bright central 
part of a galaxy that shines via the energy released 
by materials that falls toward the super massive black 
hole (SMBH). During its active phase, the AGN lumi¬ 
nosity can outshine that of the host galaxy, emitting 
an enormous amo unt of energy (lO^^-lO^® erg s“^; e.g., 
iWoo fc Urrvll200^ over a wide range of wavelengths from 
gamma-ray to radio. The spectrum of an AGN provides 
us with important information concerning the structure 
and the physics around the SMBH that powers the AGN 
activity. High-energy photons, such as X-ray and UV, 
probe the accretion disk around the SMBH that gives 
off hot thermal emission. Optical photons trace the ac¬ 
cretion disk or the broad and narrow line regions farther 
away from the SMBH. Meanwhile, infrared (IR) observa¬ 
tions can trace the dust-heated radiation from the dusty 
torus or the star forming regions. Finally, jets can be 
probed in the radio. 
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In recent years, the near-infrared (NIR) wavelength re¬ 
gion has become an important window for understanding 
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Additionally, the NIR wavelength region includes in¬ 
teresting emission/absorption line features such as the 
hydrogen Paschen and Brackett lines, the polycyclic aro¬ 
matic hydrocarbon (PAH) emission feature at 3.3/xm, 
and the absorption lines from H2O ice (3.1/xm), CO2 
(4.26^m), and CO (4.67^m) gas. These line features 
can be used to understand the BLR size, the BH mass 
(Mbh), the star formation, and the materials surround¬ 
ing the nuclear region as outlined below. 

Recently, NIR Paschen lines hav e been put for¬ 
ward as a usefu l Mrh estimator ()Kim et al.l 120101 
iLaudt et all [200! IWl . since they are much less af¬ 
fected by dust extinction than the common Mrh esti¬ 
mato r s using UV/optical lines (iVestergaard fc PetersonI 
I2OO6I: iMcLure fc Dunlop 12004 iGreene fc Hoi 1200*517 A 
similar study can be extended to the Brackett lines which 
would be even less affected by the dust than the Paschen 
lines. In the case of E{B — U) = 2 mag (a value corre¬ 
sponding to the typical color excess of dusty, red AGNs; 
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e.g.. iGlikman et ^120071 : fUrrutia et alJl2009li . the Bra, 
Br/3, Pa, and P/3 line fluxes would be suppressed by fac¬ 
tors of 1.31, 1.62, 2.16, and 3.97, respectively, assum¬ 
ing t he galactic extinction law with Ry = 3.07 (iMcCalll 
I2004f) . In other words, the Bra line is three times less 
affected by the dust extinction than the P/3 line in this 
example. 

The 3.3 /rm PAH feature h as been suggested as a possi- 
ble star formation indicator (|Tokunaga et al.iri991ll , and 
has been detected in many local AGNs. A number of 
studies suggest that there may exist a correlation be¬ 
tween the starburst (traced by the 3.3 /xm PAH feature) 


and the AGN activity (traced, e.g., by 5100 A continuum 
luminosity, hereafter, A 5100 )) although the co rrelation 
may break down at ve r y high IR luminosities (lOi et al. 
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2012 

: lYamada et al.ll2013(). Finallv. molecular lines such 


as H 2 O ice, CO 2 , and GO arise from molecular clouds 
or the diffused in terstellar medium in the line of sight 
toward the AGN (iSpoon et al.ll2004) . and they are use¬ 
ful tracers for the reservoir of molecular gas in the host 
galaxy and in the vicinity of an AGN. NIR hydrogen line 
ratios are also suggested as a use ful means to und erstand 
the physical conditions of BLR (iRuff et al.ll2012ll . 

Studies of Brackett lines and PAH features have a 
special importance for understanding the evolutionary 
stages of galaxies, especially the dusty AGN phase. 
Dusty, red AGNs are AGNs whose optical light is ob- 
scured by the foreground dust fsee iKim et al.l l2?)Tfll for 
the demography of various types of red AGNs). They 
are suspected to be an intermediate population in the 
evolution of galaxies, situated between the dust en¬ 
shrouded star forming phase (such as ultra luminous in¬ 
frared galaxy, ULIRG) and the luminous AGN phase. 
During the red AGN stage, SMBHs are expected to grow 
rapidly p roviding feedback to the star forma tion of host 
galaxies (iHopkins et al.ll2005t iLi et al.ll2007tl . but heavy 
dust extinction in host galaxies prohibits observational 
studies of the activities occurring in red, dusty AGNs. 
Brackett lines and PAH features can possibly serve as 
NIR spectral diagnostics that can explore the AGN and 
star formation activity of dusty, red AGNs due to the 
low extinction of the light in NIR. 

While studies of NIR spectra at A < 2.5/xm are com¬ 
mon, studies based on the 2.5-5.0/xm spectra of nearby 
AGNs have been scarce, mainly due to the observa¬ 
tional limitation caused either by atmospheric absorp¬ 
tion or strong thermal background at A > 2.5 /xm on the 
ground. In order to better understand the 2.5-5.0/xm 
spectra of AGNs, we conducted a spectroscopic study 
of 83 near by AGNs using the AKARI Infrared Ca m- 
era (IRG) (jMurakami et al.l[2007t lOnaka et all 120071) as 
a part of the AKARI mission program. Quasar S pectro - 
scopic Observation with NIR Grism lQSONG: llml201(Tl l. 
The AKARI IRG covers the wavelength range 2.5-5.0 /xm 
with R ^ 120 at 3.6 /xm and offers high sensitivity at 
these wavelengths thanks to the cold and transparent 
environment in space. The scientific goals of our study 
are to (1) provide a spectral atlas of nearby, bright AGNs 
that have been studied extensively in previous works; ( 2 ) 
establish the Mbh estimator using Brackett lines so that 
we can estimate Mbh reliably by reducing the dust ex¬ 
tinction effect (the Brackett-based Mbh estimator will be 


an improvement over Paschen line-based estimators); (3) 
understand the correlation between the hot dust emis¬ 
sion and other AGN properties; and (4) study the 3.3 /xm 
PAH feature in AGNs to understand the connection be¬ 
tween AGN and star-formation activities. In this paper, 
we will describe the sample characteristics, present the 
2.5-5.0 /xm AKARI spectra of 83 AGNs, derive the Mbh 
estimator using the hydrogen Brackett lines of 10 AGNs 
where the Brackett lines are well detected, and investi¬ 
gate the correlation between the hot dust emission and 
AGN properties. The investigation of the PAH feature 
will appear in a future work. 

2. THE SAMPLE AND OBSERVATION 
2 . 1 . The Sample 

When constructing the sample, we impose the follow¬ 
ing two criteria. In order to place one of the redshifted 
Brackett lines (Br/3 and Bra) in the IRG spectroscopic 
wavelength range of 2.5-5.0 /xm, our sample is limited to 
AGNs at z < 0.5. In addition, we selected objects that 
are brighter than 14 mag in the K band to gain S/N > 
10 for the continuum detection. 

Using the above criteria, we constructed a sample 
of 108 nearby, bright AGNs that are culled from the 
three different AGN source catalogs below. The first 
source of the sample is a catalog of 35 bright and 
nearby type-1 AGNs with Mbh v alues measured from 
the reverberation rna.pping method (iPeterson et al.ll2004 
iDennev et al.l 120101 : I Grier et all 120081) . The reverbera- 
tion mapping method provides the reference Mbh val¬ 
ues, which are used as the basis for deriving Mbh es¬ 
timators using spectra taken at a single-epoch. Among 
the reverberation mapping sample, 31 AGNs were ob¬ 
served before AKARI stopped its operation in 2011. 
Note that two objects among these 31 AGNs, 3G 273 
and NGG 7469, were observed through another AKARI 
mission program (AGNUL; PI: Takao Nakagawa). The 
AGNs that are selected from the reverberation mapping 
sample have a wide range of Mbh values, 1 O®'^- 1 O® ®M 0 , 
and they are at z = 0.002 - 0.29. Th e second source is the 


PG QSO catalog (|Greei] 
optically very bright an 
in previous works (iKasi 

et alJll986n. since PG QSQs are 
d have been extensively studied 
pi et al.l 120001: IVestereaardl 120021: 

IVestereaard & Peterson 

120061). We find 69 PG QSQs at 


z = 0.03 - 0.48 to satisfy our selection criteria, excluding 
those that are already included in the reverberation map¬ 
ping sample. Among them, 48 were observed by AKARI. 
For their Mbh values, we use the values derived from 
the H,5 line as listed in fVestergaard &: PetersonI (I2006D 
where the Mbh values span the range of 10®-^-10 ®-^Mq. 
Finally, we also supplement the sample by four nearby 
bright quasars from Seoul N ational University Quasar 
Surve y in Qptical (SNUQSQ; llm et al.l l2007t iLee et al] 
I 2 OO 8 I) . SNUQSQ is a survey of bright quasars missed by 
previous quasar survey at high and low galactic latitude 
(b < I 20° I), and it provides additional bright QSQs to 
the sample. In total, 108 bright, low-redshift QSQs are 
selected, of which 83 were observed by AKARI 
Figures 1 and 2 show the basic characteristics of the 
sample. AGNs in our sample are located at a red- 
shift range of z = 0.002 - 0.48, and span over a wide 
range of Mbh (10®'^-10®'^Mq), bolometric luminosity 
(Lboi; erg s“^) and accretion rate (-2.35< 
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Figure 1. Top: redshift distribution of AGNs in our sample. The 
red, blue, and green histograms are for PG QSOs, reverberation 
mapped AGNs, and SNUQSOs, respectively. Bottom: redshift vs. 
A'-band magnitude of our AGNs where the meaning of the colors 
is identical to the top panel. The A'-band magnitudes include the 
host galaxy light. Note that PG QSOs are about two magnitude 
fainter than the reverberation mapped AGNs in apparent magni¬ 
tude since they are in general at higher redshifts. 


log(Eddington ratio) <0.43). The Lboi values of the PG 
QSOs and reverberation-mapped AGNs are estimated 
from Lflinn using a b olometric correction factor of 10.3 
(iRichards et al.l 1200611 . and their L,»sinn values are take n 
from iBentz et al.T l 20091) : [Vestergaard fc Pet erson| ( 2006 b 
For SNUQSOs, t he Ha luminosity from llm et al.l ( 2007 1: 
iLee et a hl (I2008D is c onver ted to L 5100 using Equation 
01 of iGreene fc (l2005h . and the converted L 5100 is 
translated to Lboi with a correction factor of 10.3. All of 
the objects observed by AKARI in our sample and their 
basic properties are presented in Table 1. We note that 
the Mbh values in Table 1 a re calculated wit h a recently 
updated virial factor of 5.1 (jWoo et al.ll2013[1 . 


2.2. The Observation 

The 2.5-5.0/im NIR spec tra were obtained with the 
IRG infrared spectrograph (jOnaka et al.l l2007t) on the 
AKA RI satellite during pha se 3 of the AKARI oper¬ 
ation ([Murakami et al.| l2007f) . i.e., our observation was 
carried out during 2008-2011 after the helium cooling 
was over. We used the NG grism mode which provides 
a spectral resolution of R = 120 A/3.6 /rm. The spectral 
resolution corresponds to the FWHM velocity resolution 
of 2500 km s“^ or tr = 1062 km s“^ at 3.6/rm which is 
adequate for studying the broad emission lines of quasars 
and AGNs whose typical FWHM velocity widths are 
greater than 2000 km s“^. Also, the wavelength cov¬ 
erage of 2.5-5.0^m enables us to sample Br/3, Bra, and 
PAH at low redshift. 

For most of the targets, three-pointing observation was 
performed where one-pointing observation is made of 
seven to nine frames with 44 s of exposure each. For 
a few sources, 1, 2, or 5 pointings were assigned (total 


integration time of 396 to 1848 seconds) due to the ob¬ 
servational constraints of the telescope. Each target was 
placed at the center of a square-shaped slit window (l.'O 
X l.'O size) designed to avoid confusing the spectrum with 
neighboring objects. Table 1 summarizes the observation 
of each target. 

2.3. Data Reduction 

The data reduction was performed using the stan¬ 
dard data reduction software package “IRG Spectroscopy 
Toolkit for Phase 3 ^ ta Version 20110301 (2011/03/01)” 
(|Ohvama et al.ll2?i(y^ . but with an additional procedure 
for cosmic ray/hot pixel removal to improve quality of 
the final spectrum. First, we used the standard IRC 
pipeline on each pointing data separately, up to the point 
before the extraction of the one-dimensional spectrum. 
Then, we removed hot or bad pixels in the stacked two- 
dime nsional spectrum using L.A.Cosmic (jvan DokkumI 
12001 in spectroscopic mode, since the AKARI phase 
3 images typically come with many hot pixels and we 
wanted to remove effectively the hot pixels. We set the 
hot or bad pixel detection limit to 4.5tT for removal of 
the hot pixels. Then, a one-dimensional spectrum of 
each pointing for each object was extracted using an ex¬ 
traction aperture width of 5 pixels, for which aperture 
correction is performed using the pipeline and assum¬ 
ing a point source PSF to correct for the missing flux 
due to the adoption of the finite aperture size. The ex¬ 
traction aperture width gives the optimal S/N since the 
FWHM of PSFs in the IRC data is about 5 pixels (7'.'3). 
Note that the AGN spectrum can be contaminated by 
the host galaxy light, since the square-shaped slit win¬ 
dow (l.'O) does not completely block the light from the 
extended regions of a host galaxy. Hence, we are ef¬ 
fectively sampling the light from the AGN-|-host galaxy 
within an observational rectangular aperture of 7'.'3 x 
Ho, which corresponds to physical scales of 1.49 kpc x 
12.3 kpc at z = 0.01 and 32.5 kpc x 267 kpc at z = 0.3 
where the long axis runs along the dispersion direction. 
In both cases, the aperture contains a significant portion 
of host galaxies, although we expect that the measured 
Brackett lines are dominated by emission from the nu¬ 
clear region, as we describe in Section 4.1. Finally, the 
extracted one-dimensional spectra from all of the point¬ 
ings of a source were combined. In this process, any 
remaining outlier pixels (hot or bad pixel features) in 
the one-dimensional spectrum were removed by a clip¬ 
ping with 4.5 ct where the median value and the cr of a 
data point were determined including two adjacent data 
points. The wavelength and flux calibration was per¬ 
formed in the IRG pipeline and the expected accuracies 
of the wavelengths an d fluxes are about 0.01 /rm and 10% 
(jObvama et al.l[2007h . respectively. 

To demonstrate the accuracy of the photometric cal¬ 
ibration, we compare the fluxes of PG QSOs from the 
AKARI spectra to WISE W1 and Spitzer IRAG Chi 
magnitudes. To simplify the comparison, we restrict 
the sample for the comparison to point-like sources (i.e., 
unresolved) in all of the AKARI, WISE, and Spitzer 
images. Thus, we restrict the sample for the pho¬ 
tometry comparison to those in the PG QSO sample. 
For this comparison, we used the WISE W1 magni¬ 
tudes measured with a profile-fitting method provided 
by the WISE all-sky source catalog. However, for the 
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Figure 2. Left: Mbh values vs. values of AGNs in our sample. The meaning of the symbols are identical to Figure 1. Our sample 

covers a wide range in Mbh ' Mq) and Lbol erg s“^). Right: Mbh values vs. Eddington ratios of our AGNs. Our 

AGN sample covers an extensive range of BH accretion rates (-2.35 < log(Eddington ratio) < 0.43). Note that PG QSOs and SNUQSO 
AGNs are intrinsically more luminous and have higher accretion rates than the reverberation mapped AGNs. 
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Figure 3. Left: comparison of the WISE W1 or Spitzer IRAC Chi magnitudes vs. the WISE W1 and Spitzer IRAC Chi magnitudes 
derived from the AKARI spectra. The blue circles and the red squares indicate the WISE W1 and the Spitzer IRAC Chi fluxes, 
respectively. The black solid line is for the case where the AKARI ^ and WISE or Spitzer fluxes are identical. Right: Distribution of the 
ratios of the WISE or Spitzer fluxes to the AKARI fluxes of PG QSOs. The blue and the red histograms indicate the ratios for WISE 
and the Spitzer fluxes, respectively. 
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Spitzer IRAC Chi magnitudes, we measured the mag¬ 
nitudes using IRAC Chi PBCD imaging data that we 
obtained from the Spitzer archive. In order to mea¬ 
sure the magnitudes from the IRAC Chi PBCD imag¬ 
ing data, we used the IRAF daophot task with an 
aperture radius of SxFWHM. Among the PC QSOs, 
five objects (PC 0007-f 106, PC 0050-f 124, PC 1100-f772, 
PG1501-fl06, and PG1704-f608) have IRAC Chi 
PBCD imaging data. Excluding PC 0050-1-124 which is 
too bright and saturated at the center of the IRAC image, 
we used four objects from the Spitzer archive for the com¬ 
parison. Finally, using AKARI spectra, we calculated 
the corresponding WISE W1 and Spitzer Chi magni¬ 
tudes by applying appropriate filter transmission curves. 
Figure 3 shows the comparison of the fluxes from AKARI 
with those from WISE and Spitzer. The AKARI fluxes 
are on average 10% larger than the WISE fluxes, while 
the Spitzer Chi fluxes are in good agreement with the 
AKARI flux values. Although the reason for the sys¬ 
tematic offset in the flux is not clear, we find the 10% 
difference to be acceptable consid ering the general flu x 
calibration uncertainty of AKARI (lOhvama et al.l[2007f) . 

The flux errors of the spectra from the pipeline appear 
somewhat larger than the fluctuation of the data points 
(e.g.. Figure 7). To check the flux error values from the 
pipeline, we used independent flux measurements derived 
from several pointing observations (Section 2.2), and we 
checked the determined uncertainty of the flux at each 
wavelength by adopting the rms value of the flux mea¬ 
surements from each pointing data. We find that the 
error derived this way agrees with the AKARI pipeline 
result, suggesting that the AKARI pipeline error mea¬ 
surement reflects unknown fluctuations in the continuum 
level determination or sensitivity fluctuations of IRC. 

2.4. Construction of Composite Spectrum 

In an attempt to increase the S/Ns of the AGNs 
without Brackett line detections, we created a com¬ 
posite spectrum of the 48 PC QSOs. The PC QSOs 
are less affected by host galaxy contamination than the 
reverberation-mapped sample due to bright luminosity 
(Lboi: erg s“^; Figure 2). 

To construct the composite spectrum, we stacked, 
deredshifted, and flux-calibrated spectra of PC QSOs in 
logarithmic scale. Considering PG QSOs have different 
luminosity ranges, a geometric mean was used to con¬ 
struct the composite spectrum. We built a wavelength 
grid with a pixel resolution twice the spectral resolution 
of IRC (R = 120 A/3.6 fim). In each wavelength grid, the 
weighted mean flux of each object was stacked for the 
geometric mean flux. 

Figure 4 shows the composite spectrum of PG QSOs in 
A versus Fi,. The location of key lines are also displayed 
in the figure. The composite spectrum is available in 
the supplemental material accompanying the article. Al¬ 
though we could not detect Brackett lines in individual 
PG QSO spectra, some of the NIR Hydrogen emission 
lines are now visible in the composite spectrum. 

3. AKARI NIR SPECTRA 

We show the AKARI spectra from our program in Fig¬ 
ure 5. The machine-readable form of the spectra in Fig¬ 
ure 5 are available. Several interesting lines are also in¬ 
dicated in Figure 5, such as the Hydrogen Brackett lines. 



2.0 2.5 3.0 3.5 4.0 4.5 

A [pm] (Rest Frame) 


Figure 4. Top: composite spectrum of the 48 PG QSOs in ar¬ 
bitrary Fi/ units. The composite spectrum indicates Pa, Br/3, and 
PAH lines. Bottom: the number of spectra contributing for the 
composite spectrum. 


NGC4151 



Figure 6. NIR spectrum (2.5-5.0/^m) of NGC4151. The 2.5- 
5.0 /i.m wavelength window includes Bra and Br/3 lines and various 
molecular lines such as H 2 O, [Si IX], CO 2 , [Mg IV], and CO. The 
blue and the red solid lines indicate the continuum fit to the re¬ 
gions around the Br/3 and Bra lines, respectively. The dashed lines 
indicate the wavelength range used for the continuum fitting. 

3.3 //m PAH, and other molecular/atomic lines. Bra and 
Br,5 lines are identified in 9 and 7 objects at S/N > 3, 
respectively, and we mark the location of these lines in 
Figure 5 when they are identified. Potential detections 
are marked with a “?” sign. On the other hand, the 
Paschen lines are identified in 5 objects. All of the Brack¬ 
ett line detections come from the reverberation mapping 
sample whose members are brighter than the QSOs in 
the other two samples. The rather low detection rate of 
the Brackett lines reflect the reduced sensitivity and the 
increased number of detector defects during the warm 
mission phase of AKARI. 

4. BRACKETT LINES 

Using the spectra we presented in the previous sec¬ 
tion, we measured the line fluxes and widths of Brackett 
lines and used them to estimate the Mbh values and the 



























6 


Kim et al. 


physical conditions of the gas in BLR. In this section, 
we will describe the measurements of the line fluxes and 
widths, present the Brackett-line based Mbh estimators, 
and discuss the BLR physical condition inferred from the 
line luminosity ratios of the NIR Hydrogen lines. 

4.1. Brackett Line Luminosity and Width 

In this subsection, we describe how the Brackett lines 
are fitted to derive line fluxes and widths. The fitting of 
the lines starts with the identification of the line in the 
deredshifted spectra. As mentioned already in the pre¬ 
vious section, we find Bra and Br/3 at S/N > 3 in 9 and 
7 AGNs, respectively. After the line identification, the 
continuum around the Brackett lines were determined 
by fitting a linear function to both sides of each Brackett 
line. Although a NIR spectrum of a quasar is dominated 
by the black body radiation from the hot part of the dust 
torus, the linear fit of the continuum should provide a 
reasonable approxi mation of the local continuum aroun d 
each emission line (|Landt et al.l l2008t iKim et ^ 1201011 . 
However, determining the best wavelength region for the 
continuum is not a simple task, because many molec¬ 
ular absorption or atomic emission lines exist around 
the Brackett lines. Such lines are H 2 O (2.53, 2.66, 2.73, 
2.75-3.55 /im), [Si IX] (3.93 /im), CO 2 (4.26/rm) and CO 
(4.67 /im). Hence, we chose the continuum regions that 
are least affected by such lines. The wavelength regions 
adopted for the continuum determinations are listed in 
Table 2. 

After continuum subtraction, we fit the Brackett lines 
and molecular/atomic lines simultaneously. Fitting of 
the molecular/atomic lines near the Brackett lines is nec¬ 
essary, since the fitting of the Brackett lines can be in¬ 
fluenced by the presence of molecular/atomic lines such 
as the H 2 O absorption lines at 2.53 and 2.66/rm (for 
Br/3), and the [Si IX] and CO 2 gas line at 3.93 and 4.1- 
4.4/tm (for Bra). Figure 6 shows a 2.5-5.0/tm spectrum 
of NGC4151, which includes Bra, Br/3, II 2 O, [Si IX], 
and CO 2 lines. The Blue and red solid lines are the con¬ 
tinuum levels used for the Br/3 and Bra emission line 
measurements. The region bracketed by the dashed lines 
indicates the wavelength range used for the continuum 
fitting of each Brackett line. 

We used a single Gaussian function to fit each line con¬ 
sidering the low spectral resolution of the AKARI spec¬ 
troscopy, setting the central wavelength, the line width, 
and the line flux to be free parameters. A n interactive 
data language (IDL) procedure, MPFITEXPR ([Markwardtl 
l2nMi . was used to fit the lines. The fit provides the line- 
of-sight velocity dispersion, a, and fluxes of the Brackett 
lines. Note that the narrow emission line does not sig¬ 
nificantly affect the line flux and a measurements. As 
a test, we set a model spectrum with a broad (a: 850 
- 2500 km s“^) and a narrow line component (a: 170 
km s“^) at the AKARI spectral resolution. We adopt 
the flux ratios of the narrow to broad component flux 
to be 0.II5 and 0.227 for Br/3 and Bra, respectively. 
These ratios are derived by first taki ng a median Ha nar - 
row/broad component flux ratio in iLandt et al.l (l2008f l 
for five objects (MrkllO, Mrk79, NGG7469, NGC4593, 
and NGG4I5I), and then converting the Ha flux to 
the Brackett flux based on th e computation from the 
GLOUDY code (version 10.00; iFerland et al.l for 


the BLR condition (a = -1.0, n = 10® cm and U = 
10“^ ®; IKim et all 1201011 and the narrow line condition 
(n = cm“® and U — 10“®; lGoxll2000ll . We fitted 
the degraded model emission line with a single Gaussian 
function. The measured flux and a values are 108% (104 
- 115%) and 95% (90 - 98%) compared to the original 
broad component properties, which would affect Mbh by 
a small amount, ^0.03 dex. It is also confirmed that the 
narrow emission lines in the NIR do not affect the line 
flux and a measurements as much (< 10 %) in a relatively 
low-resolution spectra ('e.g.. lKim et al.l[2010ll . 

For further analysis in the following subsections, we ex¬ 
clude objects with substantial a errors (>60% of the mea¬ 
sured a value) or measured Brackett line FWHM (2.35 a) 
less than the spectral resolution (R = 120A/3.6/rm). 
The excluded cases are MrkllO (Br/3) and Ark 120 
(Bra). Finally, eight Bra and seven Br/3 lines are se¬ 
lected and used for further investigation of the line prop¬ 
erties. All of these are from the reverberation mapping 
sample, since the PG QSOs are generally too faint for a 
reliable fitting of the lines (e.g., Figure 1). 

Note that the selected Brackett lines show a small sys¬ 
tematic shift in the wavelength of 5\ = (Ameasured — Aiab) 
~ —0.014/im ((dAura) = —0.015/im and (dAur/?) = 
—0.013/im). The Bra and Br/3 lines are shifted by a 
similar amount in wavelength (dAnra — ^Aur/? = — 0.001 ± 
0.005), suggesting that the wavelength shift applies to the 
entire spectral range. The shift is a bit larger than the 
expected wavelength offset er ror in the AKARI IRC of 
0.01 /im ([Ohvama et al.l[2007T l. The origin of the system¬ 
atic offset is unclear, although it seems to be instrumen¬ 
tal. A potential source of the offset is molecular/atomic 
absorption affecting the peak wavelength of the emission 
lines. However, if this were true, then we would expect a 
stronger offset in Bra than in Br/3 since Br/3 is generally 
much less affected by molecular absorption features than 
Bra, but this was not the case. 

The formal fitting errors are found to be about 27% 
(between 7% and 49%) in tr and 28% (between 6 % and 
46%) in flux, due to low S/N and the necessity of simul¬ 
taneously fitting the other lines in the vicinity. Another 
source of error comes from the determination of the con¬ 
tinuum. We varied the wavelength regions for the con¬ 
tinuum fit from the values noted in Table 2 to see how 
strongly the changes in the continuum level affect the 
line fitting results. We find that the variation in the con¬ 
tinuum level changes the derived a and flux values by 
13% and 8 % respectively. Overall, combining these two 
errors, we obtain the error in a and flux to be 33% and 
30%, which corresponds to 0.30 dex in Mbh- 

Energy sources of the Brackett lines could be BH ac¬ 
tivity or star formation activity. Considering the large 
slit aperture of the IRC and the extraction width we 
adopted, the measured Brackett line fluxes may include 
the star formation contribution from the host galaxy. 
To determine the dominant energy source of the Brack¬ 
ett lines, we compared the expected Anra from star 
formation to the me asured Lurry The relationships 
T 3 . 3 PAH/T 1 R ~ 10“^ (|Mouri et al. Il99fll; llmanishil 1200211 
and LBra/-^^iR 10 “^ (jlmanishi et al.l I 2010 H are ex¬ 
pected for starburst galaxies, suggesting that the lumi¬ 
nosity of the Bra line from star formation is 10% of 
the 3.3 /im PAH line luminosity. Among our seven Bra- 
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NGC 7469 NGC 7469 



X [|j,m] (Rest Frame) 


Figure 7. Gaussian fits of the Brackett lines of our sample. The black solid lines and bars indicate the observed spectra in rest frame 
and the associated errors. The blue and red lines indicate the fits for Br/3 and Bra, respectively. The dashed lines show the fits to the 
molecular/atomic lines around Brackett lines, such as H 2 O, [Si IX], and CO 2 lines. The solid lines indicate the Brackett line component of 
the best-fit spectra. 
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Figure 8. Comparison of the line luminosities of the Brackett lines vs. the H/3 line. The dotted line indicates the best-fit correlation 
between the luminosities of the H/3 and the Brackett lines. 
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Figure 9. Comparison of the cr widths of the Brackett lines vs. those of the H/3 line. The dashed line indicates the case where the cr 
value of the H/3 and the Brackett lines are identical. The weak correlation for Bra is likely to be due to low S/N and the CO 2 absorption 
feature right next to the Bra line. The CO 2 absorption is the most prominently seen in the spectrum of NGC4151, which is represented 
by the red filled circle. 
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detected AGNs, the 3.3 /xm PAH emission line is detected 
in only two objects (NGC4051 and NGG 7469) using the 
extraction width adopted in this study (Kim et al. in 
prep). The expected star formation contributions of the 
Bra line based on the 3.3 /xm PAH line are 8% and 46% 
for NGG 4051 and NGC 7469, respectively. This assumes 
the above 10% flux ratio of Tera to Ts.spah, but the re¬ 
lation contains a large scatter of a factor of two to three. 
Therefore, we conclude that the contribution of star for¬ 
mation to the Bracket emission lines is negligible (no or 
weak detection of PAHs) for most of the AGNs in our 
sample. NGG 7469 is an exceptional case with a very 
strong 3.3/xm PAH. Even in this case, the star forma¬ 
tion contribution is less than half for Bra, but we cau¬ 
tion possible over-estimation of the line flux (134%) and 
under-estimation of the line width (80%) for NGC 7469. 

The measured values of the tr and luminosities of the 
Brackett lines and the corresponding H/3 properties are 
presented in Table 2. Note that the a values are corrected 
for the instrumental resolution by subtracting the instru¬ 
mental resolution in quadrature (tr^ = — uf^gt) from 

the observed values. The derivation of the line luminosi¬ 
ties assumes a standard ACDM cosmol ogy of Ho=70 k m 
s“^ Mpc“^, r2m=0.3 and flA=0.7 ie.g.. Ilm et all[19^ 


4.2. Correlation between Balmer and Brackett Lines 

In this section, we will investigate how the properties 
of the Brackett lines correlate with those of the Balmer 
lines. If the two lines arise from BLRs similar to each 
other, then we expect their line properties to correlate 
well. This would justify th e use of the Bra ckett lines 
as an Mbh estimator (e.g., IKim et all I2010H . just like 
the Balm er lines that are common ly used as Mbh es¬ 
timators ([Greene fc Hoi [20051 1200711 . For this compar¬ 
ison, we take the average of the multi-epoch a values 
and the luminosities of the H/3 line from previous studies 
(iKaspi et al.ll2000l: iMarziani et al.l 120031: iPeterson et al.l 

12004 IVestergaard fc Petersonl I2006t iLandt et al.l 1200811 . 
When there is a large discrepancy in the a values (> 
xl.5) of H/3 and Ha, we take the average of the two. 
Such a case occurred for only one object (PG 1411-1-442). 
Note that we could not find the H/3 line luminosity 
data for NGC 4051, s o we only adopted a a value from 
IPeterson et al.l (1200411 for NGC 4051. 

Figures 8 and 9 show the correlations between the 
Brackett lines and the H/3 line in luminosity and a. In 
Figure 8, we compare the luminosities of the Brackett 
lines with those of the H/3 line. The Pearson correlation 
coefficients are 0.991 and 0.947 for Br/3 and Bra, respec¬ 
tively. By performing a linear bisector fit to these points, 
we find that the two quantities correlate with each other 
as 


+ (0.85±0.05)tog( io/°;^^-. ). 

( 1 ) 


and 

+ (0.79 ±0.07) log 

( 2 ) 

In Figure 9, the dashed line indicates a perfect correla¬ 
tion between the a values of the Brackett and H/3 lines. 
The figures show a reasonable correlation between the 
Br/3 and H/3 line widths, but a weak correlation between 
the Bra and H/3 line widths. The Pearson correlation 
coefficients are 0.735 and 0.273 for Br/3 and Bra, respec¬ 
tively. The weak correlation in a for Bra is likely due to 
low S/N and the difficulty in fitting some of the Bra lines. 
The zodiacal background increases rapidly at A > 4 /xm, 
causing the background noise to increase for Bra. Fur¬ 
thermore, the existence of the CO 2 absorption feature 
right next to Bra makes it challenging to provide a good 
fit. CO 2 absorption is the most prominent in the spec¬ 
trum of NGC 4151, which also seems to be a significant 
outlier in Figure 9.The exclusion of the NGC 4151 point 
improves the a correlation for Bra, with the correlation 
coefficient of 0.468. 

We also compared the Brackett line luminosities to the 
BLR radii. For the radii of BLR, we adopted values 
from the time l ags in the va r iabilit y between the H/3 and 
the L 5100 from iBentz et al.l (|2013| 1. Figure 10 shows the 
correlation between the Brackett line luminosities and 
the BLR radii. The best-fit relations are 

+ (0.62±0.05)log( i^„^«^-, ), 

(3) 

and 

+ ( 0 . 66 ± 0 . 06 )log( i^,„^°;°^-, ). 

(4) 

4.3. BH Mass Estimators with Brackett Lines 

In this section, we present Mbh estimators that are 
based on the Br/3 and Bra lines. The idea here is that 
the radius and the velocity of the BLR can be inferred 
from the luminosity and a of the Brackett lines, similar to 
the Mbh estimators based on the Balmer of the Paschen 
lines. Mathematically, we need to find three unknown 
parameters a, b, and c in the following equation: 

log(M) = a-I-61og(L)-I-clog(cr). (5) 

For simplicity, we fix the exponent of the velocity term, 
c to 2, as expected in the virial theorem. The exponent 
of the luminosity term, 6, is expected to be 0.5 theo¬ 
retically and emp i rically, i.e., Rbt.r ~ (iDibail 119771 : 
iKaspi et al.l l200(]l : iGreene fc Hd 120051 : IKim et al.l 1201011 . 
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Figure 10. Comparison of the line luminosities of the Brackett 
lines vs. the BLR sizes. For the BLR size, we use the time lag in the 
variability between the H/3 and -L 5100 from IBentz et anil2013ll . The 
solid line indicates the best-fit correlation between the luminosities 
of the Brackett lines and the BLR sizes. 


However, we will treat 6 as a free parameter. In order to 
deriv e the Mrh estimator s, we used the MPFITEXY rou¬ 
tine (iWilliams et al.l[2010ll to fit a linear relationship in 
the logarithmic scale as in Equation (5). MPFITEXY re¬ 
turns the fitted parameters with the uncertainties of the 
parameters and an intrinsic scatter in the fitted relation. 
The routine has been tested against several other sim- 
i lar routines, and has shown to return reliable outputs 
(jPark et al.ll2012ll . As a result, we obtain the following 
relations: 


/ O-Br/3 

V10^ km s“^ / 

( 6 ) 

/ ^Bra 

V 10^ km s“i / 

In Figure 11, we plot the Mbh from the reverberation 
mapping method versus the Mbh from the Brackett-line- 
based estimators, i.e., the values from Equations (6) and 
(7). The rms scatters of the data points (including mea¬ 
surement errors) with respect to the best-fit correlations 
of Equations (6) and (7) are 0.229 dex (Br/3) and 0.407 
dex (Bra). The intrinsic scatter of Equation (7) is 0.284 
dex, but the intrinsic scatter of Equation (6) is not de¬ 
termined since the fit did not require intrinsic scatter. 
The value of the intrinsic scatter is expected to be small 
(< 0.2 dex). 

Since the cr values are rather poorly determined from 
the Brackett lines but we expect that they correlate 
closely with the a values of the H/3, we also derived 
the Mbh estimators using L^r:a,p but adopting cth^ as 
a proxy for a B rn.a- The cth/? values are taken from 
iPeterson et al.l (1200411 . Since dust extinction in a host 
galaxy can reduce its line luminosity significantly but 
not as much for the line width, this is a plausible path 
to estimate Mbh of red AGNs whose can be deter¬ 
mined to a good accuracy, but its H/3 luminosity is quite 
uncertain due to an unknown amount of dust extinction. 

In this process, we also fixed the exponent of the veloc¬ 
ity term, c, to 2. The resultant Mbh estimators using 


-^BH _ j^q6.61±0.23 

Mq 

and 

-^BH _ ^q6.68±0.20 

Mq 


LbtP 


104O gj.g g-1 


^Bra 

1040 gj.g g-1 


0.68±0.13 


0.66±0.21 


Asra,/? and a^p are given below and also in Figure 12 as 


Mbh 

1^ 

and 


10 


6.44±0.17 


LbtP 


1040 gj.g g 


-1 


0.72±0.10 


( \ 

V10^ km / 

( 8 ) 


A/T / T \ 0.48±0.17 ^ 

-'^BH _ , „6.76±0.17 ( ABra \ / CTBp 

Mq \10^°ergs“^/ VlO^kms"^ 

(9) 

The raw rms scatters (including measurement er¬ 
rors) with respect to a perfect correlation between the 
reverberation-mapping-based Mbh and Mbh for Equa¬ 
tions (8) and (9) are 0.203 dex (Br,5) and 0.364 dex 
(Bra), while the intrinsic scatters are 0.152 dex (Br/3) 
and 0.310 dex (Bra). However, the intrinsic scatters are 
not better than the values for the purely Brackett-based 
Mbh estimators in Equations (6) and (7). 

Since the a values of Brackett lines, especially for Bra, 
show a weak correlation, we test whether or not the cr 
values of Brackett lines are adding useful information 
to the Mbh estimators. For this test, we derive Mbh 
estimators using only Brackett line luminosities: 


-^BH _ ^q6.92±0.26 
Mq 


L 


Br/3 


and 


-^BH _ ^q6.79±0.15 
Mq 


104O gj.g g-1 

Abfo 

1040 gj-g g-1 


0.73±0.15 


0.76±0.14 


( 10 ) 


( 11 ) 


’ Figure 13 shows the Mbh values from the reverbera¬ 
tion mapping method versus the Mbh values from the 
Mbh estimators using only LbtP or Lnra, i-C., the values 
from Equations (10) and (11). We find that the intrinsic 
scatters in these relations are 0.314 dex (Br/3), and 0.251 
•dex (Bra). Therefore, for Br/3, the intrinsic scatter in the 
Mbh estimators becomes smaller as we move from the 
luminosity-only estimator to the full Br,5 luminosity-Fcr 
estimators, lending support to the usefulness of Equation 
(6). On other hand, the same cannot be said for the Bra- 
based estimators. This suggests that Bra cr measurement 
is more challenging with the current data, and the eval¬ 
uation of the Bra-based Mbh estimators is needed with 
data sets with better sensitivity. Considering that the 
Eddington ratio of AGNs used for Equations (10) and 
(11) are between 0.015 and 0.225, Equations (10) and 
(11) could be useful when (i) the Brackett line luminosi¬ 
ties are available but not for reliable cr values and (ii) the 
Eddington ratio of the sample is within a range of 0.015 
to 0.225. 


4.4. Line Luminosity Ratios 

We investigate the physical conditions of BLR by 
comparing the observed line luminosity ratios of the 
Balmer through Brackett lines to the line luminosity ra- 
tios computed from the CLOUDY code (version 10.00; 
iFerland et al.l[lM^ . The line luminosity ratios are com¬ 
puted by varying three physical parameters - the shape 
of the ionizing continuum (a = -1.0, -1.5), the ioniza¬ 
tion parameter {U = 10“®'^ - 10°'®), and the hydrogen 
density (n = 10® - 10^^ cm ^). The hydrogen number 



































11 





Figure 11. Comparison of the Brackett-line-based Mbh vs. Mbh derived from the reverberation mapping method. The Brackett-line- 
based Mbh is derived by using the L and a of the Brackett lines as Equations (6) and (7). The dashed line indicates a line where Mbh from 
Brackett lines and the reverberation mapping methods are identical. In each panel, we indicate the intrinsic and the raw (no correction of 
the measurement errors) scatters of the points with respect to the dashed line. 





Figure 12. Comparison of Mbh from Equations (8) and (9) vs. Mbh derived from the reverberation mapping method. The Equations 
(8) and (9) are established by using the Z/Bro:,^ but adopting cth/S- The meaning of the line is identical to Figure 11. 



Figure 13. Comparison of Mbh from Equations (10) and (11) vs. Mbh derived from the reverberation mapping method. The Equations 
(10) and (11) are established by using only Z/Bra,/?- The meanings of the line is identical to Figure 11. 


densities are c ommonly taken as n 10^ - 10^° cm~^ for 
type-l AGNs (iDavidson fc NetzeHll979l : iRees et ^119891 : 

iFerland &: Perssonlll989ll . but there are studies suggest¬ 
ing a high hydrogen number density (n > 10^^cm“^) 
with a low ion ization parameter of 17 = 10“^ (e-g-, 
IRuff et a1.l[2nT2ll for BLR. 

Figure 14 shows the observed line luminosity ratios 
(with respect to H/3) as a function of wavelength for 
seven AGNs, for which we have at least one Brackett 
to Balmer line ratio and an additional line ratio of ei¬ 
ther Paschen or Brackett lines. The line luminosities 
are corrected fo r galactic ext inction with the York Ex¬ 
tinction S olver (iM cCalll 12 0041). adoptin g E{B — V) val¬ 
ues from (|Schlaflv fc Finkbeineiil201lll and Ry = 3.07. 
No attempt is made to correct for internal dust ex¬ 


tincti on. Note t hat th e Paschen line ratios are taken 
from iKim et ^ (j201fl[ l. The mean line ratios of the 
seven AGNs are iera/TBr/s = 1-37 ± 0.31, Lb^p/Lup = 
0.073 ± 0.010, and TBra/Tn/? = 0.081 ± 0.014. 

In order to find the most likely physical condition for 
each object, we searched for a parameter set that min¬ 
imizes x^, which is a function of line ratios such as 
Lpp/Lbp, Lpry/Lua , LBra/Lp a , and Aurry/Tufl follow¬ 
ing the approach by IRuff et al.l (|2012[1 : 


2 _ X ^ (77observed,z Rmodel^i) 

^ / j 2 

1 — 1 ^ 


( 12 ) 


where N is the number of line luminosity ratios, and 
two types of Ri are the line luminosity ratios from either 
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Figure 14. Line luminosity ratios of Paschen and Brackett with respect to H/3. The crosses and the circles represent AGNs in the 
reverberation mapped sample with the Brackett line lumin osity measuremen ts. When available, we also plot the line luminosity ratios for 
Pa (squares) and P/3 (triangles) with respect to H/3 from fKim et al.l H2010IL The gray lines represent the expected line ratios from the 
CLOUDY code, assuming various parameters with different combinations of a = —1.0 —1.5, n = cm~^, and U = 

In each panel, we indicate the most likely set of the parameters, with the CLOUDY model predictions for such a parameter set in the red 
dashed line. 


Table 3 
BLR Conditions 


Object Name 

a 

n 

[cm“®] 

u 

Mrk 79 

-1.5 

10“ 

10-1.5 

Mrk 110 

-1.5 

lOio 

10-0.5 

NGC4151 

-1.5 

10® 

10-5.5 

3C273 

-1.0 

lO^i 

10-1.5 

NGC4593 

-1.0 

lO^i 

10-2.5 

PG 1411-1-442 

-1.5 

10® 

10-1.5 

NGC 7469 

-1.5 

loll 

10-5.5 


observation (when with a subscript, “observed”) or the 
CLOUDY model (with a subscript, “model”), and ai is 
the error in the measured line ratio. The best-fit param¬ 
eter sets are indicated for each object in Figure 14, and 
are summarized in Table 3. The best-fit line luminosity 
ratios are also plotted as a red dashed line in Figure 14. 

The CLOUDY models can reproduce the observed line 
luminosity ratios, under commonly cited BLR physical 
conditions of n = 10® - 10^^ cm“^ and U = 10“^ ® - 
10® ®, suggesting that most of the objects with Brackett 


line detections do not require a more exotic condition 
(n > 10^® ^ < 10”^)) such as that found 

in I Ruff et al.l (I2012D . One exception is NGC7469, for 
which we find a best-fit parameter set of n = 10^^ cm“® 
and U = 10“®'®. This is the case where the Brackett 
line luminosity may be affected by the star formation 
component (possibly by as much as 50%), and thus the 
interpretation of the line luminosity ratio result cannot 
be taken too seriously for this object. 


5. DUST COMPONENT 

In this section, we present hot ( > 1000 K; e.g., 
iBarvaini^ Il987t iG likman e t al.l 1200 6) and warm 
200 K; e.g.. lNetzer et al.ll2007tlDeo et al.l[20?)^ dust tem¬ 
peratures and luminosities of the PG QSOs and the 
reverberation-mapped AGNs, and we investigate the cor¬ 
relation between hot dust luminosity (Lhd) and Lboi- 
For this purpose, we restrict our analysis to luminous 
AGNs with Lboi > 10^®'^ergs“^ for which the host 
galaxy light cont ribution is less than 10% at 0.51 /iin 
(jShen et al.ll201ll l. and even smaller at NIR. 

To study the AGN dust components, we assembled a 
multi-wavelength data set from the optical to MIR by ob¬ 
taining photometry of the objects from SPS S, 2MASS, 
WISE (ITS and 1F4), and ISO (7.3//m: iFIaas et ^ 





























13 


NGC 4593 


C2. 



1.00 


0.10 


0.01 


a=-1.0 log(n)=11. log{U)=-2.5 - 



Bra 


12 3 4 


NGC 7469 


02 . 



1.00 


0.10 


0.01 


a=-1.5 log(n)=11. log{U)=-5.5 - 



Brp Bra 


12 3 4 

X [|im] 


Figure 14. Continued 

l2003f) . These multi-wavelength photometric data are 
combined with the AKARI spectra, and we fit them with 
a SED model that is composed of a single power law and 
a double black body radiation component as 

F, = Coiy" + CiS,(Thd) + C2B,(Twd), (13) 

where a is the slope of the underlying power-law contin¬ 
uum which dominates the optical light, Fv is the Planck 
function as a function of Thd and T-wd are the hot 
and warm dust temperatures, and Cq, Ci, and C 2 are 
the normalization constants of each component. The fit¬ 
ting result gives us the luminosities and the temperatures 
of the hot/warm dust components. 

We note three technical points that we took into ac¬ 
count during the fit. First, we excluded SDSS photom¬ 
etry between 0.6 and 0.7 ^m in the rest frame to avoid 
the inclusion of the Ha line which generally has a large 
equivalent width and produces a bump in the SED that 
deviates from a power law which is assumed as the con¬ 
tinuum shape in the optical. Second, we increased the 
errors of the SDSS data points to represent the deviation 
of the optical continuum shape from a power law. Here, 
the deviation is caused by the intrinsic fluctuations in the 
continuum spectra due to Fe emission complex, emission 
lines other than Ha, and the host galaxy light. With¬ 
out increasing the error, a fit to the SDSS photometry 
points with a simple power-law function ends up with a 
rather large reduced values due to very small original 
SDSS photometry errors. From a composite spectrum of 
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a quasar (|Glikman et al.l 1^00611 placed at 0 < z < 0.4, 
we find the intrinsic fluctuation of the optical spectrum 
with respect to a simple featureless power law to be 
Gm = 0.035 mag, and this value is added in quadrature to 
the original SDSS photometry error in each band. Third, 
we calculate the reduced chi-square, xt of the fit, and the 
hot and warm dust properties are taken seriously only if 
xl < 5. One object (PG1322-1-659) is excluded by the 
X^ limit for which we had difficulty fitting the shape of 
the AKARI spectrum (the last object in Figure 15). 

Through the procedure above, we obtain the tempera¬ 
tures and luminosities of the hot and warm dust com¬ 
ponents of 19 AGNs (2 reverberation mapped AGNs 
and 17 PG QSOs). The fitting results are presented in 
Figure 15, showing that both the hot and warm dust 
components exist and contribute to the NIR and MIR 
spectra of AGNs. However, four objects (PG 1114-1-445, 
PG1116-b215, PG1543-f489, and PG1545-b210) show 
an excess of 7.3 ^m ISO photometry. This excess can 
hint at the existence of an intermediate-temperature dust 
component or a continuous dust temperature. 

The SED fit provides formal fitting errors of 0.8% in 
Thd and 3.1% in Lhd, where we obtain a small error 
in Thd due to the strong constraints provided by the 
AKARI spectra. The errors in Table 4 indicate the 
formal fitting errors. Since the changes in the assump¬ 
tions that went into the fitting model could influence 
the result more than the formal fitting errors, we ex- 
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Figure 15. Photometric and spectroscopic SED. The SEDs are shown in rest-frame. Triangle, circle, square, and star indicate photometric 
data point from SDSS, 2MASS, WISE (VE3 and 1^4), and ISO (7.3/j-m), respectively. The black solid line denotes AKARI spectrum. The 
dashed lines are fitted results. The blue, red, and yellow dashed lines mean power-law, hot, and warm dust components, respectively, and 
the black dashed line is the sum of all components. The measured hot and warm dust temperatures are indicated in the lower right-hand 
corner. The dust properties are taken only if < 5, and the x5 each sample is indicated at top of each figure. The PG 1322-1-659 is 
excluded by the xj limit, and the SED of the PG 1322-1-659 is indicated in last panel of this figure. 


Table 4 

Dust Gomponents Fitting Results 


Objects 

Thd 

[K] 

Two 

[K] 

log(I/HD) 
[erg s-l] 

log(LwD) 

[erg s-l[ 

CFhd 

a. 

PG 1211+143 

864+5 

191+2 

45.27+0.01 

45.38+0.02 

0.28 

-1.55+0.02 

PG 1307+085 

854+13 

184+3 

45.27+0.03 

45.45+0.03 

0.18 

-1.49+0.02 

PG 0043+039 

1068+14 

230+9 

46.22+0.02 

46.04+0.07 

0.23 

-1.83+0.03 

PG 0947+396 

1036+2 

213+3 

45.78+0.00 

45.68+0.03 

0.60 

-1.68+0.01 

PG 1048+342 

1032+20 

213+5 

45.06+0.03 

45.12+0.04 

0.25 

-1.55+0.04 

PG 1049-005 

1068+6 

222+2 

46.48+0.01 

46.65+0.02 

0.43 

-1.68+0.03 

PG 1100+772 

1062+9 

215+3 

46.20+0.02 

46.14+0.03 

0.29 

-1.89+0.04 

PG 1103-006 

1144+17 

250+5 

46.33+0.03 

46.37+0.03 

0.35 

-1.95+0.05 

PG 1114+445 

917+2 

195+2 

45.56+0.00 

45.73+0.02 

0.51 

-1.50+0.01 

PG 1116+215 

1078+1 

216+3 

46.17+0.00 

45.84+0.03 

0.37 

-2.00+0.02 

PG 1121+422 

1078+15 

236+9 

45.55+0.02 

45.11+0.06 

0.21 

-1.95+0.04 

PG 1259+593 

1326+1 

315+8 

46.81+0.00 

46.32+0.04 

0.41 

-2.15+0.01 

PG 1402+261 

1106+3 

219+2 

45.95+0.00 

45.79+0.02 

0.64 

-2.23+0.01 

PG 1543+489 

1231+2 

239+2 

46.58+0.00 

46.81+0.02 

0.83 

-1.98+0.01 

PG 1545+210 

1151+8 

209+3 

46.03+0.00 

45.82+0.03 

0.29 

-2.17+0.04 

PG 1704+608 

1145+3 

191+2 

46.63+0.01 

46.67+0.03 

0.36 

-1.55+0.01 

PG 2112+059 

1172+7 

263+4 

46.91+0.00 

46.82+0.03 

0.52 

-1.75+0.03 

PG 2233+134 

1111+7 

224+3 

46.08+0.01 

46.26+0.02 

0.40 

-2.25+0.05 

PG 2308+098 

1142+10 

177+17 

46.29+0.02 

46.08+0.15 

0.19 

-1.50+0.03 



































15 


100.0 


10.0 


1.0 


0.1 


100 


E 10 


A SDSS 
• 2MASS 
■ WISE 


PG 1048+342 

Xv= 0.66 




Thd=1032K 

Tw=213K 


1 10 

PG 1100+772 



1 


10 


1000 i-— 

: A SDSS 
• 2MASS 
■ WISE 
100 ^ ★ ISO 


10 : 


PG 1114+445 

Xv= 2.97 


Tho=917K 
!„□= 195K 


1000E 


100 : 


10 : 


A SDSS 
• 2MASS 
■ WISE 


PG 1049-005 

Xv= 1-94 


Thd=1068K 
T,„n= 222K 


1 10 

PG 1103-006 


100F 


10 : 


A SDSS 
• 2MASS 
■ WISE 


Xv= 107 



Thd—1144K 
T„b=250K X- 


1 


10 


1000 i—— 

: A SDSS 
• 2MASS 
■ WISE 
100 ^ * ISO 


10 : 


PG 1116+215 

Xv= 3.38 




Thd=1078K 
T„b= 21 6K 


1 10 1 10 
Wavelength [|j.m] Wavelength [jim] 


Figure 15. Continued 

















































16 


Kim et al. 


100E 


E 10 : 


A SDSS 
• 2MASS 
■ WISE 


PG 1121+422 

Xv= 1-37 


-f' 


PG 1259+593 


Thd=1078K 
!„□= 236K 


1 10 
PG 1402+261 


1000 


100 


10 


A SDSS 
• 2MASS 
■ WISE 
★ ISO 


Xv= 1-52 




10 


PG 1545+210 


100 


10 



1 10 
PG 1543+489 


lOOOf 


100 ; 


A SDSS 
• 2MASS 
■ WISE 
★ ISO 


Xv= 3.00 




-^^ 


10 

“ vT 

- !- -^- 1 

Thd=1106K 

'5 



Tho=1231K 

T™=219K 




239K 


10 


PG 1704+608 


IOOf 


E 10 : 


: - - - 1 

- 1 


1UUU 


' 


A SDSS ■A= 1.53 



: A SDSS 2.04 


• 2MASS 



• 2MASS 



■ WISE 




■ WISE 

_ 


★ ISO 

J - 


100 

- ★ISO 


_ 






A--'*'/ 





10 




' ~ T ^ ! 

Tho=1151K 

\ " 

^' 7 

Thd=1145K 


* 1 - - - 

!„□= 209K 

V 



T„o=191K 





1 




1 

10 



1 

10 


Wavelength [|j.m; 




Wavelength [gm; 




Figure 15. Continued 


plored how the derived quantities could be affected by 
changing the following conditions in the fit as an attempt 
to estimate more realistic errors: (1) fixing the power- 
law slope of the continuum, a, in Equation (13) to a 
value that deviates by five times the formal error from 
the ht; (2) the reduction of the fluxes of the AKARI 
data by 10%, since AKARI fluxes could be overesti¬ 
mated by that amount (Figure 3), and (3) the subtrac¬ 
tion of the host galaxy contribution. The mean absolute 
differences in the fitted quantities from the change in 
the power law slope are 3.3% (Thd; between 0.6% and 
8.7%) and 11% (Thd; between 2.6% and 36%), and from 
the flux of AKARI data are 2.7% (between <1% and 
6.5%) in Thd and 15% (between 11% and 22%) in Thd- 
The reduction in the AKARI flux leads to the reduc¬ 
tion in Thd and Thd- We note that the change of the 
power-law slope and the AKARI flux value keeps the xl 
value to less than five, except for one object. We esti¬ 


mated the contribution of the host galaxy light in T nmn 
(|Peterson et al.l[200il: IVestergaard fc Peters(^ 12006(1 us¬ 

ing Equation (li of iShen et al.ll2011l and subtracted the 
host galaxy light assuming an elliptical galaxy SED shape 
(|Bruzual fc Charlotl I2003f) with the conditions of solar 
metallicity, r = 0.1 Gyr, t = 10 Gyr, and a Chabrier 
IMF. We examined how the fitting result changes with 
the subtraction of the host galaxy component, and find 
mean differences of 1.5% (between <1% and 7.5%) in 
Thd and 7.5% (between <1% and 22%) in Thd, where 
Thd decreases and Thd increases in most cases. With 
these considerations, the error of Thd is realistically 
^4%, while the error in Thd could be as much as ^14%. 
For the warm dust component, we have only two pho¬ 
tometric data points for most AGNs and a single point 
in some AGNs; we caution that the htted result have 
larger errors and can be affected more easily by the ht- 
ting method than the hot dust component. 
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Figure 15. Continued 


The left panel of Figure 16 shows Lboi versus Thd- 
The Lboi values are estimated fro m Lsinn using a bolo - 
metric correction factor of 10.3 (iRichards et al. [200®, 
and the Tmnn values are tak e n from iBentz et al.l 1)20091 1 
and iVestergaard fc Petersoiil (|200® . The Lboi and the 
Lhd values have a positive correlation as 


Thd 
1800 K 


2^q-0.231±0.003 


Lboi 


1046 gj-g g 


-1 


0.061±0.005 

(14) 


and the Pearson correlation coefficient is 0.690. We find 
that the mean Thd and Twd are 1083 K and 221K, 
respectively, which i s lower than the commonly cite d 
Thd of 1500 K (e.g., IBarvainiill987t lElvis et HI 1 199411 . 
and somewhat lowe r than the value qu oted in a recent 
study 11260 K. from lGlikman et afll20ci® . The difference 
does not seem to arise from technical parts of the fitting 
method. Whe n we measure the Thd using a compos¬ 
ite spectrum of lGlikman et al.) (|200® with our method, 
we find the result to be identical to that of the previous 
study (i.e., 1260 K). Also, we estimate Thd values o f our 
sample using the methods of lGlikman et al.l (1200® . and 
the measured Thd values used in this way are consistent 
with our results. 

The right panel of Figure 16 shows the hot dust lumi¬ 
nosities versus the Lboi values for the 19 AGNs. For the 
bolometric luminosity, we use L 5100 from our SED fit- 
ting multiplied by a b olometric correction factor of 10.3 
()Richards et al.l l200® . In the figure, we also indicate 
lines with three different values for the covering factor 
of the hot dust component (CFhd)- Here, the cover¬ 
ing factor is given by the ratio of the luminosity from 


the hot dust compon ent to Lboi (CFHD=LHD/Lboi; e.g., 
iMaiolino et al.ll200^ . We find that Lhd correlates well 
with Lboi, and their mean covering factor is found to be 
0.38. The measured me an covering factor is similar to 
several previous results (ISanders et al.lll989li lElvis et HI 
II 994 iRoseboom et al.l 1201® . However, we note that 
there is no significant downturn in CFhd with Lboi, in 
contrast to previous results which reported that CFhd 
has an anti-correlatio n with Lboi at the high luminosity 
end of 10^^^-^ erg s~^ (iWang e t al.l 1200,^ IMaioli no et all 
20071: iTreister et 120081 : iMor fc Trakhtenbrod 120111 : 
Mor fc Netzei]l2012f) . This discrepancy is quite possibly 
due to the low number statistics in our data, consider¬ 
ing that the anti-correlation found by other groups has 
a large scatter. 


6. SUMMARY 

Using the AKARI IRC in grism mode, we obtained 
2.5-5.0^m spectra of 83 nearby (0.002< z <0.48), bright 
{K <14 mag) type-I AGNs. The Brackett lines are iden¬ 
tified in II AGNs in the sample, allowing us to derive 
Mbh estimators that are based on Bra and Br/3 line 
properties. We find the Br/3-based Mbh estimator us¬ 
ing the measured line width and luminosity so as to be 
useful with an intrinsic scatter of ~0.I dex or less with 
respect to the Mbh values from the reverberation map¬ 
ping method. On the other hand, the Bra-based esti¬ 
mator shows a larger scatter and no strong dependence 
on (TBra, demonstrating the difficulty of CBra measure¬ 
ments and the need for future data to refine the esti¬ 
mator. However, both Mbh estimators can potentially 
be applied to derive Mbh values of dusty AGNs due 
to the low extinction in the Brackett line wavelengths. 
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Figure 16. Left: Lboi vs. Thd of 19 AGNs. The solid line indicates the best-fit correlation. Right: comparison between L\^o\ 
Lhd- The solid, dotted, and dashed lines denote CFhd = 0-5, and 0.1, respectively, and the measured mean CFjjd of fho 19 AGNs is 

0.38. 


Table 5 

Spectrum of Mrk 335 


A 

[mH 

fx 

[erg s ^ cm ^ /i.m 

f\ Uncertainty 
[erg cm“^ 

2.543 

3.500E-11 

4.653E-13 

2.552 

3.363E-11 

4.595E-13 

2.562 

3.207E-11 

3.482E-13 

2.572 

3.133E-11 

2.153E-13 

2.581 

3.106E-11 

3.536E-13 

2.591 

3.118E-11 

5.891E-13 

2.601 

3.128E-11 

2.988E-13 


Note. — This table display only a part of spec¬ 
trum of Mrk 335. The entire spectra of 83 AGNs 
and the composite spectrum of 48 PG QSOs are 
available as a tar file in the electronic version of 
Astrophysical Journal Supplement Series. 

The continuum shape is well traced for all of the objects, 
and we derived the temperatures and luminosities of the 
hot and warm dust components of 19 luminous AGNs 
(Lboi > lO"^^'^ erg s“'^) for which host galaxy contami¬ 
nation is expected to be small (< 10%). The use of the 
AKARI spectra provides the advantage of tracing the 
NIR continuum densely, thus improving the dust tem¬ 
perature measurements. Our results show the hot dust 
temperature to be about HOOK, which is lower than the 
commonly cited value of 1500 K but broadly consistent 
with most recent hot dust temperature measurements. 
The covering factor of the hot dust component is also de¬ 
rived with a mean CFhd of 0.38, consistent with previous 
measurements. The AKARI spectral atlas and the tab¬ 
ulated spectral measurements can be downloaded from 
this article, revealing the rarely studied spectral range of 
2.5-5.0/xm for low-redshift AGNs. 
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7. APPENDIX 

We provide reduced 2.5-5.0^m spectra of 83 AGNs 
and the composite spectrum of 48 PG QSOs. Table 5 is 
an example spectrum of Mrk 335. The full version of the 
reduced spectra is available in machine readable table 
form. 
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Figure 5. Spectra and images of AGNs in the QSONG sample. The left panels show de-redshifted 2.5—5.0/im spectra and the spectral 
regions affected by remaining hot or bad pixels are masked out. The ordinate is in units of erg s“^ cm~^ . Several hydrogen 

lines (Pa, Br/3, and Bra) and molecular lines (PAH, CO 2 , and CO) are marked on the spectra. When the line detection is not obvious, 
the line is marked with a question mark. The right panels indicate object images from the digitized sky survey (DSS). The white dashed 
line box and the cyan solid bar indicate I'O x I'O slit window and 7'/3 width of extraction aperture, respectively. The arrow at the top 
left denotes north. 
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Table 1 
Object List 


Objects 

a 

(J2000.0) 

S 

(J2000.0) 

Redshift 

z 

K 

(mag) 

log(A4'BH) 

(Me) 

Exposure Time 
(=5) 

Observed Dates 

Mrk 335 

00 06 19.5 

+20 12 10 

0.025 

10.59 

7.11^ 

1188 

2008-06-30 

PG 0026+129 

00 29 13.6 

+13 16 03 

0.142 

12.11 

8.56=" 

1100 

2008-07-02, 03 

PG 0052+251 

00 54 52.1 

+25 25 38 

0.155 

12.23 

8.53=" 

1144 

2008-07-13, 14 

Fairall 9 

01 23 45.7 

-58 48 20 

0.047 

11.11 

8.37=" 

1188 

2008-12-01, 02 

Mrk 590 

02 14 33.5 

-00 46 00 

0.026 

10.57 

7.64=" 

660 

2009-07-24 

3G120 

04 33 11.0 

+05 21 15 

0.033 

10.57 

7.71=" 

1188 

2009-02-25 

Ark 120 

05 16 11.4 

-00 08 59 

0.032 

10.22 

8.14=" 

1056 

2008-09-09 

Mrk 79 

07 42 32.7 

+49 48 34 

0.022 

10.68 

7.68=" 

1188 

2008-10-10 

PG 0804+761 

08 10 58.6 

+76 02 42 

0.100 

10.91 

8.80=" 

1188 

2008-10-04 

PG 0844+349 

08 47 42.4 

+34 45 04 

0.064 

12.01 

7.93=" 

1188 

2009-04-24 

Mrk 110 

09 25 12.8 

+52 17 10 

0.035 

12.25 

7.36=" 

1056 

2008-10-29, 2009-04-25 

NGC3516 

11 06 47.4 

+72 34 06 

0.008 

9.61 

7.46'" 

1188 

2008-10-25, 26 

NGC 3783 

11 39 01.7 

-37 44 18 

0.009 

10.01 

7.44=" 

1188 

2008-07-03 

NGC4051 

12 03 09.6 

+44 31 52 

0.002 

10.01 

6.20'" 

748 

2009-05-29, 30 

NGC 4151 

12 10 32.5 

+39 24 20 

0.003 

8.51 

7.09=" 

1056 

2008-06-03, 04 

PG 1211+143 

12 14 17.7 

+14 03 12 

0.080 

11.29 

8.13=" 

1144 

2009-12-18, 19 

3C 273 

12 29 06.7 

+02 03 08 

0.158 

9.97 

8.91=" 

1848 

2008-06-26, 27, 12-27 

PG 1229+204 

12 32 03.6 

+20 09 29 

0.063 

12.14 

7.83=" 

1188 

2009-06-20 

NGC 4593 

12 39 39.4 

-05 20 39 

0.009 

10.03 

6.69=" 

1144 

2008-07-02, 2009-01-01 

PG 1307+085 

13 09 47.0 

+08 19 48 

0.155 

12.42 

8.61=" 

1144 

2008-07-03, 04 

Mrk 279 

13 53 03.4 

+69 18 29 

0.030 

10.39 

7.51=" 

1144 

2008-11-14, 15, 16 

PG 1411+442 

14 13 48.3 

+44 00 14 

0.089 

11.50 

8.61=" 

1144 

2008-06-28, 29, 30 

NGC 5548 

14 17 59.5 

+25 08 12 

0.017 

10.13 

7.61'" 

1144 

2008-07-13, 14 

Mrk 817 

14 36 22.0 

+58 47 39 

0.031 

10.90 

7.60=" 

660 

2008-06-12 

PG 1613+658 

16 13 57.1 

+65 43 09 

0.129 

11.41 

8.41=" 

1144 

2008-06-07, 08 

PG 1617+175 

16 20 11.2 

+17 24 27 

0.112 

12.41 

8.74=" 

1188 

2008-08-21 

PG 1700+518 

17 01 24.8 

+51 49 20 

0.292 

11.98 

8.85=" 

924 

2008-08-15, 17 

3C 390.3 

18 42 08.9 

+79 46 17 

0.056 

11.71 

8.42=" 

1188 

2008-09-10 

Mrk 509 

20 44 09.7 

-10 43 24 

0.034 

10.19 

8.12=" 

1012 

2009-04-30 

PG 2130+099 

21 32 27.8 

+10 08 19 

0.062 

10.66 

7.54"= 

1100 

2009-05-20 

NGC 7469 

23 03 15.6 

+08 52 26 

0.016 

9.762 

7.05=" 

1100 

2008-06-10, 11 

PG 0003+158 

00 05 59.2 

+16 09 49 

0.450 

13.53 

9.24=' 

1144 

2008-06-30, 12-29, 30 

PG 0007+106 

00 10 31.0 

+10 58 30 

0.089 

11.78 

8.69=' 

1188 

2009-12-28 

PG 0043+039 

00 45 47.3 

+04 10 24 

0.385 

13.59 

9.09=' 

1012 

2010-01-02 

PG 0049+171 

00 51 54.8 

+17 25 58 

0.064 

13.34 

8.31=' 

1188 

2010-01-09 

PG 0050+124 

00 53 34.9 

+12 41 36 

0.061 

10.35 

7.41=' 

1144 

2008-07-08, 09 

PG 0838+770 

08 44 45.3 

+76 53 10 

0.131 

13.15 

8.12=' 

1188 

2008-10-06 

PG 0923+129 

09 26 03.3 

+12 44 04 

0.029 

11.30 

8.56=' 

396 

2009-11-12 

PG 0934+013 

09 37 01.0 

+01 05 43 

0.050 

13.50 

7.01=' 

792 

2009-05-17 

PG 0947+396 

09 50 48.4 

+39 26 51 

0.205 

12.76 

8.64=' 

1188 

2009-11-07 

PG 1011-040 

10 14 20.7 

-04 18 40 

0.058 

12.65 

7.28=' 

1144 

2009-11-28 

PG 1012+008 

10 14 54.9 

+00 33 37 

0.186 

13.08 

8.21=' 

1100 

2009-11-27 

PG 1022+519 

10 25 31.3 

+51 40 35 

0.044 

12.71 

7.11=' 

396 

2009-11-07 

PG 1048-090 

10 51 29.9 

-09 18 10 

0.344 

13.95 

9.17=' 

1188 

2009-06-08 

PG 1048+342 

10 51 43.9 

+33 59 27 

0.167 

13.60 

8.33=' 

704 

2009-11-22 

PG 1049-005 

10 51 51.4 

-00 51 18 

0.359 

12.84 

9.15=' 

1144 

2009-06-04 

PG 1100+772 

11 04 13.7 

+76 58 58 

0.311 

13.05 

9.24=' 

1188 

2008-10-18 

PG 1103-006 

11 06 31.8 

-00 52 52 

0.423 

13.85 

9.29=' 

792 

2009-06-08 

PG 1114+445 

11 17 06.4 

+44 13 33 

0.143 

12.33 

8.56=' 

1144 

2009-11-21, 22 

PG 1115+407 

11 18 30.3 

+40 25 54 

0.154 

12.76 

7.63=' 

1144 

2009-11-24 

PG 1116+215 

11 19 08.7 

+21 19 18 

0.176 

11.54 

8.48=' 

1056 

2009-12-03 

PG 1121+422 

11 24 39.2 

+42 01 45 

0.225 

13.26 

8.00"' 

396 

2009-11-26 

PG 1202+281 

12 04 42.1 

+27 54 12 

0.165 

12.87 

8.58=' 

396 

2009-12-11 

PG 1216+069 

12 19 20.9 

+06 38 39 

0.331 

13.30 

9.16=' 

616 

2009-12-24 

PG 1244+026 

12 46 35.2 

+02 22 09 

0.048 

13.19 

6.49=' 

1100 

2008-06-30, 12-31 

PG 1259+593 

13 01 12.9 

+59 02 07 

0.477 

13.07 

8.88^ 

1144 

2009-11-27, 28 

PG 1302-102 

13 05 33.0 

-10 33 19 

0.278 

12.85 

8.84=' 

1100 

2009-07-10, 11 

PG 1322+659 

13 23 49.5 

+65 41 48 

0.168 

12.85 

8.25^ 

1056 

2009-11-21, 22 

PG 1351+236 

13 54 06.4 

+23 25 49 

0.055 

13.01 

8.53=' 

792 

2008-07-09 

PG 1402+261 

14 05 16.2 

+25 55 34 

0.164 

12.16 

7.91=' 

1144 

2008-07-09 

PG 1415+451 

14 17 00.7 

+44 56 06 

0.113 

12.23 

7.98=' 

1144 

2008-12-29, 2009-06-28 

PG 1416-129 

14 19 03.8 

-13 10 44 

0.129 

13.34 

9.01=' 

1100 

2010-01-27 

PG 1425+267 

14 27 35.6 

+26 32 15 

0.364 

13.57 

9.70=' 

1188 

2008-07-16 

PG 1427+480 

14 29 43.1 

+47 47 26 

0.220 

13.48 

8.05=' 

1188 

2009-06-29, 07-01 

PG 1448+273 

14 51 08.7 

+27 09 27 

0.065 

12.05 

6.94=' 

1188 

2009-07-22, 23 

PG 1501+106 

15 04 01.2 

+10 26 16 

0.036 

11.47 

8.49=' 

1188 

2009-08-02 

PG 1512+370 

15 14 43.0 

+36 50 50 

0.370 

13.69 

9.34=' 

1100 

2008-07-22 

PG 1519+226 

15 21 14.2 

+22 27 44 

0.137 

12.47 

7.91=' 

1188 

2009-08-02, 03 

PG 1534+580 

15 35 52.3 

+57 54 09 

0.029 

12.18 

8.17=' 

1100 

2008-06-29, 12-31, 2009-01-01 
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Table 1 — Continued 


Objects 

OL 

(J2000.0) 

5 

(J2000.0) 

Redshift 

z 

K 

(mag) 

log(A4BH) 

(Mq) 

Exposure Time 

(s) 

Observed Dates 

PG 1535+547 

15 36 38.3 

+54 33 33 

0.038 

11.63 

7.16'^ 

704 

2008-07-06, 08 

PG 1543+489 

15 45 30.2 

+48 46 09 

0.399 

13.14 

7.96=^ 

1188 

2009-01-17, 20, 21 

PG 1545+210 

15 47 43.5 

+20 52 17 

0.264 

13.18 

9.28=^ 

616 

2008-08-11 

PG 1626+554 

16 27 56.1 

+55 22 32 

0.133 

12.67 

8.46^ 

1188 

2010-01-22 

PG 1704+608 

17 04 41.4 

+60 44 31 

0.371 

12.43 

9.36'^ 

1100 

2008-07-23 

PG 2112+059 

21 14 52.6 

+06 07 42 

0.466 

12.85 

8.97'^ 

1056 

2009-05-13, 14 

PG 2209+184 

22 11 53.9 

+18 41 50 

0.070 

12.55 

8.73=^ 

1144 

2009-12-03, 05 

PG 2233+134 

22 36 07.7 

+13 43 55 

0.325 

13.81 

8.00'^ 

1188 

2009-12-07, 08 

PG 2304+042 

23 07 02.9 

+04 32 57 

0.042 

12.54 

8.53=^ 

748 

2009-12-12 

PG 2308+098 

23 11 17.7 

+10 08 15 

0.433 

13.23 

9.56=^ 

1012 

2008-12-14 

SNUQSO 0644+3546 

06 44 10.8 

+35 46 44 

0.077 

11.82 

8.86® 

1188 

2008-10-01 

SNUQSO 1312+0641 

13 12 04.7 

+06 41 07 

0.239 

13.20 

9.09® 

1100 

2008-07-06, 2009-01-03 

SNUQSO 1935+5314 

19 35 21.1 

+53 14 11 

0.248 

13.14 

8.61® 

1188 

2008-11-17, 2009-05-11 

SNUQSO 2127+2719 

21 27 56.4 

+27 19 05 

0.197 

12.63 

8.74® 

396 

2009-05-26 


^ The Mbh values are estimated from reverberation mapping method iPeterson et al' 1 l 2004 n . 

^ The Mbh values are estimated from reverberation mapping method JDennev et al.|| 2010 D . 

The Mbh values are estimated from reverberation mapping method iGrier et aUf^OSl) . 

^ The Mbh values are based on FWHMh/ 3 and L5100 JVestergaard Sz Peterson||200^ . 

^ The Mbh values are based on FWHMho; and Lhch jLee et al ] [200a : IIm~ernil200^ . 
















Table 2. Brackett Parameters And References 


39 






u 

u 

CJ 

u 

u 


o 

O 

u 

o 

u 



1 


o 

CO 


CS 

CO 

CO 

Tp 

O 

CO 

CM 

to 


K 

CO 

f-H 


00 

00 


iO 



b- 

CM 

iC 






N- 

00 

o 

N- 

lO 

05 

N- 

to 

O 

tP 



s 










CM 


































1 

CO 



ei 

o 


Xi 


o 

d 

o 

«9 

Si 



tc 

u 


rH 

Tp 

o 


CO 


CS 

o 

t- 

to 

CM 



o 

00 


q 

CO 

1 

C^l 

CM 

CM 

N- 

q 



0) 

C'l 






d 


d 

f-H 

d 

d 



















o 

















o 



iO 

05 

O 

iO 


00 


X 


Ch 



CO 



CO 

cs 

CO 

CO 



CO 



Q 

O 

a. 


TP 

1 

1 

O iO 

23 

1 

2 

3 

2 


hJ 



cs 



CO 


CN| 

q 


i-H 

q 






TP 



TP 


TP 

TP 


TP 

TP 

TP 








00 


t- 

00 


iO 

CO 

b- 


Ch 



00 



00 

q 

00 

00 


b- 

q 

X 


Q 

u 

'S' 


CO 



CO 

CO 

CO 

CO 


co 

CO 

CO 



<i 

1 

1 

Jj 05 

c^ .4 

1 


<i 

c4 


m 

‘ * 


00 



oO 

00 

00 

00 


b- 

05 

X 





CO 



CO 

CO 

cd 

CO 


cd 

CO 

CO 

a 







00 

OI 

b- 

o 


o 



u 


1 

CO 





TP 

tP 

o 



CO 

CO 


CQ 



TP 





d 

d 


d 

TP 

CM 



bi} 

u 


-H 

1 

1 

-H 

-n 

-H 

■H 

1 

■+1 

-H 

-H 



CD 

o 





o 

CD 

b- 

CO 



o 

q 








o 

00 




ID 




o 





cd 

cd 

d 








1 


05 

CO 



rr 

1.0 

(N 

CO 

ID 

b- 

o 

CO 

Tb 

t—i 


CO 

lO 

lO 

N- 

O 

CO 

o 

•5b 


b 

CO 

S' 

-H 

1 

1 

-H 

-H 

-H 

-H 

1 

-H 

-H 

-H 



[km 


en 

b- 



CO 

05 

cs 

o 

b- 

b- 

■Cf 

b- 

N- 

00 


lO 

00 

CO 

fH 

X 

CO 

lO 

CM 

CM 






CO 

CO 




TP 

CO 


o 

X 


cn 




00 

s- 




q 

t- 


00 

q 


Q 

O 

6 

a. 



s 

(N 




d 

CM 


CM 

CM 


iO 

1 

ot 


1 

1 

1 

1 

ot 

I 

C^ 

.4 


hJ 




N. 

b; 




q 

q 


r^ 

to 






S 

ci 




(N 

CM 


CM 

CM 











00 

to 


X 

b~ 


Ch 





1.0 




q 

q 


TP 

q 


Q 

O 

'S’ 

a. 

2 

1 

S 

00 

ci 

1 

es 

1 

1 

1 

d 

ci 

CM 

1 

CO 

1 

CM 

1 

i-H 

CM 

1 

TP 


m 





lO 




lO 

iO 


TP 

iO 






c4 

<N 




d 

CM 


CM 

CM 

'XI 

y ( 


rH 








o 

05 



IO 

m 


1 

CO 



cs 

TT 




CM 


CM* 

CO 





ri 




d 

CO 


rH 

TP 



bO 

0) 

o 

S' 

1 

-H 

00 

o 

•-H 

-H 

CO 

1 

1 

1 

-H 

-H 

05 

1 

-H 

05 

-H 

to 



o 



CO 

00 




05 

d 

00 

00 



q 

00 






1.0 

b- 




o 

05 


05 

CM 






CO 

05 




o 

rb 


cs 

iO 



1 



r}< 





Tb 

b- 


CD 

N. 


b 

CO 

S' 

1 

-H -H 

1 

1 

1 

-H 

-H 

1 

-H 

-H 



g 



to 





C>l 

r-H 


(N 

to 



s 



Oi 





(N 

to 


rb 

CO 






lO 

1.0 




b- 

to 


(N 

TP 














(N 


<D 













CS 


a 

cd 

CJ 




o 

(N 

o 

iO 

+ 

(N 

lO 

o 

05 

riC 

o 

-s 

CO 

00 

b- 

CO 

U 

1.0 

o 

rb 

O 

iO 

TP 

U 

CO 

N- 

CM 

o 

CD 

05 

ID 

'5b 

u 

■5b 

+ 

t-H 

Tp 

05 

CO 

■5b 

b~ 

O 

<u 




CO 

o 

IS 


O 

O O 

CO 

o 


O 






O 








O 


o 





Oh 








CL, 



S P s 

3 3 

3.-3 

^ 

^ K c 
o > o 
u o 

^ ^ Cm 

O O 

« '2' <u 

:2 m "S 

OP _ 'm 


determination region for the Br^; (10) luminosity and (11) FWHM of H;S. For column (10) and (11), we adopted the luminosities and 
FWHMs of H/3 from literatures and make footnotes as their literatures: ^ Kaspi et al. (2000); ^ Marziani et al. (2003); Peterson 
et al. (2004); Vestergaard Peterson (2006); ® Landt et al. (2008) 









